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TtTE* 3-body Calogero problem is solved by separation of variables 
fo^-arbitrary exchange statistics. A numerical computation of 
thK4-body spectrum is also presented. The results display new 
f^ffljrres in comparison with the standard case of bosons and 
f^a^ions, for instance the energies are not linear with the in- 
teraction parameter v and Bethe ansatz as well as Haldane's 
statistics are not verified. 
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lO Introduction 

" 5 

The] Calogero model, that is the Id system of N bosons or N 
fe+jnions interacting by the 2-body potential v(y — l)x~ 2 , is of 
physical interest: it interpolates between free bosons and free 
fergjiions |l], ||, it exhibits the fractional statistics of Haldane 
|3|, and its asymptotic wavefunction is in accordance with the 
B^Sie hypothesis J|, ||. Moreover, it belongs to a large class of 
integrablc models ||. 

£Ljhe Calogero model is known to be completely integrable for 
a long time, and this fact remains true when the particles are 
in a harmonic well. Calogero has solved himself the 2-body and 
3-body problems by separation of variables. He has also shown 
that the many-body problem is integrable and he has deduced 
its energy spectrum 0. More recently, the many-body eigen- 
statcs have been constructed from an algebra of creation and 
annihilation operators ^ ^, [l(], |llj. This algebra is immedi- 
ately generalizable to the case of Boltzmann statistics where no 
exchange symmetry is imposed to the wavefunction, however the 
ground states and their energies are not all known in this case. 

In this paper, I discuss the asymptotic behaviour of the wave- 
function when two particles coincide and I apply it to solve the 
eigenvalue problem for arbitrary statistics. I analytically con- 
struct the eigenvalues and eigenfunctions for N bodies in the 
special case v = 1/2, and next for 3 bodies in the general case. 
A numerical computation of the 4-body problem is also pre- 
sented. The results obviously emcompasse the standard solution 
for bosons and fermions. However, the states of mixed statis- 
tics display a non-trivial framework and their energies are not 
linear in v. Moreover, the particle current does not necessarily 
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vanish when two particles coincide. In conclusion, I discuss the 
Calogero model for distinguishable particles with constants i/y 
,rrii depending on the labels of the particles JjJ], [l3| . 

2 The example of the two-body prob- 
lem 

2.1 The nature of the wavefunction 

Let me focus on the relative motion of two particles of coordi- 
nates X\ and X2 in a harmonic well. Then the Calogero hamil- 
tonian is 

H = -d 2 + 1 -u J 2 x 2 + »±^±, (1) 
4 x z 

with x = x\ — xi , h = 1 and m = 1 . One will assume < v < 1 
for simplicity, the limits v = and v = 1 corresponding to the 
usual harmonic oscillator problem. The hamiltonian acts on a 
domain of definition 

D H C {ip, ip G L 2 ,Hijj e L 2 ,H+ = H} (2) 

which warrants finite matrix elements, real eigenvalues and the 
orthogonality between states of different energies. There remains 
to understand the nature of the wavefunctions ip £ Dh- 

At ordinary points x ^ 0, the wavefunction and its derivative 
may be assumed continuous as usual. Around the singular point 
x = 0, one verifies that the wavefunction must have an algebraic 
behaviour as 

ip~x u or ifi-x 1 -", (3) 

otherwise its image Hip is not a square integrablc function. Note 
that the situations v < and v > 1 will not be considered in this 
paper because the wavefunction is then divergent in general. 
The hamiltonian is self-adjoint when the boundary term 

J dx (^d 2 ^ - ?M 2 ^i) ^i9^2 Or) (4) 

is continuous at any point and vanishes at infinity for all pair of 

wavefunctions i/>i,tp2- Notice that -k *0 d "4> is merely the particle 
current. It is equivalent to impose the previous condition on the 
boundary term for all pair of basic states tpi,ip2i or the same 
condition on the particle current for all wavefunctions ip. This 
condition is satisfied at ordinary points and at infinity because 
the wavefunction and its first derivative are continous at these 
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points and vanishe at infinity. At the singular point a; = 0, one 
can always set 



iP = \x\»f(x)+x\x\- v g{x) 



(5) 



and then one has 



ipi d tp2 



(1 - 2v)(f l92 - g x fc) + x{h dg 2 +gid h) 

'X 2 \x\-^( 9l d 9 2). (6) 



A2v (h a h) 



For < v < 1 and v ^ 1/2, the continuity of the particle current 
is ensured at x = if / and g are continuous up to a given phase 
shift at this point and if df and dg are finite at this point. The 
case v = 1/2 is particular because the continuity of / and g 
can then be relaxed at x = 0.[] The choice of the overall phase 
shift is essentially arbitrary. This phase shift can be removed 
by adding a constant phase depending on the ordering of the 
particles, x > or x < 0. As a result, one can demand that / 
and g are exactly continuous at x = 0. An alternative choice will 
be presented below in relation with the notion of intermediate 
statistics. 

The differential equation theory shows that / and g are entire 
series in x These functions and their derivatives are so 

simultaneously continuous at x = 0. In the limit v — 0, the 
problem is reduced to that of a pure harmonic oscillator, the 
wavefunction and its derivative are continuous as it should and 
the usual harmonic eigenstates are necessarily recovered. 

Note that both the hamiltonian and the boundary condi- 
tions are invariant under the mirror transformation ^(x^v) — > 
sign (x) ip(x, 1 — v). Introducing the odd factor sign(x), the 
mirror symmetry inverts the exchange properties of ?p, that is, 
it exchanges Bose and Fermi statistics. 

2.2 An anyon gauge 

In analogy with the 2d anyon model , it is possible to define 
a new gauge by the gauge transformation 



(7) 



with 9 — arg(x) i.e. 9 = if x > and 9 = ir if x < 0. Let me 
name it the anyon gauge. In this new gauge, the hamiltonian is 
the same but the boundary condition at x = becomes 



4/ = x v f(x)+x 1 -»g(x) 



(8) 



with / and g continuous. The last equation is well defined if 
a cut is introduced in the complex plane z — x + iy along the 
axis x < 0, and if, on this axis, x is defined by the analytic 
continuation x = lim e ^ + |a:|e 1< - 7r ~ e ^ ) . 

In the anyon gauge, the wavefunction is multivalued when 
two identical particle interchanges are successively effected and, 
thus, the wavefunction is neither symmetric nor antisymmetric 
but it has complicated exchange symmetries. However the limits 
v = and v = 1 reproduce monovalued wavefunctions: the 
gauge transformation is the identity in the case v = and it 
multiplies the wavefunction by the odd factor sign (x) in the 
case v = 1. As a result, if the wavefunction ip satisfies Bose 
statistics, the wavefunction tp' will be bosonic at v — and 
fermionic at v = 1, and vis versa if ip satisfies Fermi statistics. 

1 Note that the regular combination \x\^^' 2 f{x~) + a; | a; | 1 / 2 f(x) = 
2<5a;>o \/\x\f{%) leads to a discontinuity as 2& x >of(x). 



2.3 Eigenvalues and eigenfunctions 

The eigenstates are obtained in a standard way. Let me choose 
the regular gauge characterized by the boundary condition (|^). 
Thus, the bosonic eigenstates are 



1> n = \x\ v Ln *{-Lux')e—^ x , E n 
and the fermionic eigenstates are 



{2n+- + v)uj, (9) 



ip„ = x\x\ v LI "(^ujx 2 )e 4<- 



(2n- 



v)u, (10) 



where n is a non-negative integer. The whole basis is clearly 
orthogonal and complete in L 2 . One remarks the mirror 
symmetry between the Bose and the Fermi states. Using 
the identities H 2n {z) = (-)"2 2 "n!L^ 1/2 (z 2 ) and H 2n+1 {z) = 
(— ) n 2 2n+1 nlzLl/ 2 (z 2 ) where H n is the Hermite polynomial, the 
harmonic oscillator solution is reproduced in the limits v = 
and v = 1, except that there is a phase shift at x = in the last 
case. 

The eigenstates of the anyon gauge are merely deduced from 
the changes \x\ u — * x v and \x\~ v — > x~ v . They interpolate 
between Bose and Fermi statistics as it should. 

3 Exact results for the many-body 
problem 

3.1 The model 

The A^-body hamiltonian is 



iV 



i=l 



E 

h3< 



v(y-l) 



(ID 



with 



The wavefunction ip has to be sought in an 



irreducible representation of the permutation group Sn , usually 
bosonic or fermionic. Let me still assume < v < 1, the limits 
v = and v = 1 reproducing the standard harmonic oscillator 
problem. The domain of definition of the hamiltonian is chosen 
as for the 2-body case : the wavefunction and its derivative are 
continuous except at Xi 



Xj where 



V> = \xijVf + Xij\xij\ v g 



(12) 



with / and g locally continuous, dyf and dijg finite (dij is a 
shorted notation for dt — dj). The case v = 1/2 is special: the 
continuity of / and g is no longer necessary and the boundary 
condition at reduced to 



(13) 



with dijh finite. 



As for the 2-body case, we can define an anyon gauge by 

iw E M < e «-0 ( 14 ) 



ip = e 



where 9ij — arg(xy ). The hamiltonian is unaltered, but the 
boundary condition at Xi ^ Xi J IS changed into 



(15) 
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The standard solution of the harmonic oscillator is then exactly 
reproduced at v = and v = 1. However, at v = 1, the anyon 
gauge ( |l4| ) introduces an antisymmetric factor sign - K Xij) 
which inverts the exchange properties of the wavefunction: if 
if} belongs to an irreducible representation of Sn, then its im- 
age rp' belongs to the conjugate representation since the anti- 
symmetric factor exchanges lines (symmetrisors) and columns 
(antisymmetrisors) commuting with the Young diagramm (pro- 
jector). As a result, when v goes over [0, 1], the particles inter- 
polate between bosons and fermions, or more generally, between 
an irreducible representation and the conjugate representation 
of the permutation group Sn ■ 

The mirror symmetry reads i\)'(xi,v) — » ip'(Xi, 1 — v) in the 
anyon gauge and ip(xi, v) — > sign (n*,j< x ij)ty{ x i, 1 — v ) m tnc 
regular gauge. In the regular gauge, the mirror transformation 
inverts the exchange properties of the wavefunction introducing 
an antisymmetric factor. 

Working with an unperturbed basis, the matrix elements of 
the interaction are divergent. In analogy with the 2d anyon 
model, the non-unitary transformation 



= n 4 ^ 



v = n 



(16) 



*,j< 



(or its mirror symmetric) leads to a hamiltonian 



N , , 

ff=E(-k+^)-E^-^) ( 17 ) 



i=l 



(or respectively its mirror symmetric) whose matrix elements 
are well defined with the principal value regularization. It is 
especially adapted to the perturbative and numerical analyses. 
Albeit it is not formally hermitian, it is so between states of the 
same unperturbed energy. 

3.2 Separation of two variables 

One defines the center-of-mass coordinate X = jr Xi and a 
radial coordinate for the relative motion r 2 = ±Y,^xl. The 
remaining coordinates constitute a system of angles denoted by 
£1. In these variables, the hamiltonian reads 



H= -— d 2 x + -Noj 2 X 2 ~-8 2 - 
2N x 2 2 r 



N-2 



A 1 



2r ^ + ^ + ^ V - (18) 



where A is an angular operator. The eigenfunctions may be 
separated as i/j — l E,(X)R(r)Y(tt). Assuming that the eigenvalue 
equation AY = XY is satisfied, the normalizable solutions for the 
center of mass 



E K {X) = H K {VN^X)e 
and for the radial part 



-±NujX 2 



r v - - - - Lin (ujr 2 )e~ 
are standard. The total energy is then given by 




(19) 



(20) 



(21) 



k and n are non-negative integers. 



3.3 Raising and lowering operators 

This section is a brief account of some recent progress. The 
main idea consists to construct a covariante derivative using the 
operator which exchanges the particles i and j ||, [Io| , 



(22) 



The adjoint of Di is merely — £)$, for Pjj and are anticom- 
muting. The commutation rules are [Di,Dj] = and 



N 



[Di, Xj] = Sij(l + ^ vP ik) - vPij 



fc=l 



One next defines the hamiltonian 



1 N 



(23) 



(24) 



in terms of creation and annihilation operators, 

a? = yf^T^=, [H,at]=±uat. (25) 

These operators are mutually adjoints. 

As regards its action on bosonic wavefunctions, the hamilto- 
nian Ti. is found to be identical with the Calogero one. As a 
result, the bosonic eigenstates of the Calogero model are easily 
obtained. The Bose groundstate 



n 

i,j< 



(26) 



is annihilated by the a i 's, its energy is ^u) + N ^ ^ uiu, and 
the excitements follow from the successive actions of the af's on 
the groundstate. However, these excitements have to be sym- 
metrized under particle exchanges. To avoid this symmetriza- 
tion, it is desirable to define a set of operators preserving the ex- 
change symmetries of the wavefunction. Separating the center- 
of-mass and the relative excitements, such a set of Raising and 
lowering operators is given by 



N 

Af = y <i 



N 



i 



A"> 2 = 2> f~ ~At ) 



N 



(27) 



The set is complete if A: = 1, 2, . . . , N . For a discussion of the 
algebra associated with these operators, see [[ll] . The main com- 
mutation rules are 

[H,Af]=±kwj£. (28) 
In this way, the bosonic excitements are directly obtained as 



jv 



n^rll 



(29) 



fc=i 



and the total energies are 



En 



N 



%3< 



N 



N{N - 1) 



-L)V, 



(30) 
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where the n^'s are non-negative integers. Here, the net effect of 
the raising operators is to multiply the groundstate by a poly- 
nomial. The particle current (||) vanishes at coinciding points 
Xi = Xj because 5 = and thus there is no particle exchange. 

As regards its action on an arbitrary function, the hamilto- 
nian Ti is a non-local operator effecting among others some in- 
terchanges on the variables of the function. It is necessarily 
different from the Calogero hamiltonian H. However, H is the 
only local operator which coincides with TL acting on symmetric 
functions. Let me define as the only local operator which 
coincides with Aj^ acting on symmetric functions. The explicit 
form of can be obtained by expanding A^ in power series 
of Pij, commuting the Py's in first position and then replacing 
Pij by unity. One finds in this way 



d 



X 



At 



Hid. 



2uj 

N - 



At 



N 

- B# 



2 



2u 



N Al) 



(31) 



(32) 



(33) 



1 



2uj N 



Af), 



et caetera [||. One can directly check that A^ and A \ respec- 
tively originate the center-of-mass excitements ( |l9| ) and the ra- 
dial excitements (EOh. 



By definition, the operator A^ does not alter the exchange 
symmetries of the wavef unctions. Moreover, in a Hilbert space 
of symmetric functions, the commutation rule J2S| ) can be im- 
mediately translated into 



±kLuA^ 



(34) 



Now, this relation involves only local operators so that it is also 
satisfied independently of the exchange symmetries of the wave- 
functions. One concludes that the A^ 's are Raising and lowering 
operators for arbitrary exchange statistics. In fact, it remains 
to prove that the action of these operators on an eigenfunc- 
tion having the correct behaviour at coinciding points gives a 
new eigenfunction having also this behaviour. Since the opera- 
tors involve derivatives, it suffices to assume that the successive 
derivatives of / and g in ( fL2| ) are all continuous at Xi = Xj . Every 
eigenfunction discussed in this paper confirm this hypothesis. 

In conclusion, the ground states of the Calogero model are 
annihilated by the A^ 's and the excitements are obtained by 
the action of Ylk=i(A^) nk on these ground states. The energy 
spectrum is linear in the quantum numbers n^. According to the 
permutation theory, a complete eigenstate basis will be obtained 
if there are Nl independent ground states: one in the Bose rep- 
resentation, another in the Fermi representation, and d ground 
states in each irreducible representation of degree d of Sn- 

3.4 Eigenstates and eigenvalues for bosons 
and fermions 

The Bose eigenstates are discussed in the preceding paragraph. 
The Fermi eigenstates can be deduced from the Bose ones by 
using the mirror symmetry. 



3.5 Eigenstates and eigenvalues at v — 1/2 

At v = 1/2, the behaviour of the wavefunction when two par- 
ticles coincide is dictated by (|l^) where the continuity of h 
is no longer required. It is possible to construct an eigen- 
state ipp,n k equals to the Bose eigenstate ip nk in the sector 
x Pl > x P2 > . . . > x PN and sets to zero in the other sectors. 
The states which are non-zero on different sectors are obviously 
linearly independent. Hence one has the eigenstate basis 

N 



i> P ,n k =5 XP1>XP2> ... >XPN U(A+r* n J\xy\ e-^2> 



fe=l 



i,j< 



with the energy spectrum 



N(N + 1) 
\,n k = ^kn k iv H oj. 



(35) 
(36) 



fc=i 



The quantum numbers are the permutations p and the non- 
negative integers n^. 

Note that the energies are completely degenerated with re- 
spect to particle exchanges because the action of the permuta- 
tions on any eigenstate gives AH independent eigenstates. This 
last property implies the completness of the basis (|35|) in L 2 pro- 
vided that the subspace of the bosonic eigenstates is complete. 

Physically, at v — 1/2, the particle current always vanishes 
when two particles coincide, for all statistics. There is no ex- 
change of particles and no change of the ordering of the parti- 
cles as in classical mechanics In this situation, the coupling 
strength v{y — 1) is at its minimal value —1/4. Indeed, for a 
smaller strength, v is a complex number, the energies are also 
complex and this indicates that the particle system collapses. 

4 General solution for three particles 

4.1 System of coordinates 

Introducing the following Jacobi coordinates, 



X 

y 



1 



(xx +x 2 + x 3 ), 



1 

i 

71 



{xi 



X2) 



(37) 



(xi + x 2 - 2x 3 ), 



the relative coordinates x, y are expressed as functions of the 
standard polar coordinates r, (j> |Q] . The radial coordinate also 
satisfies the alternative definition r 2 = \ V . , x 2 -. This system 
of coordinates is related to a discrete Fourier transformation of 
the coordinates on the line, 



re 



(j xi +]X 2 + x 3 ) 



(38) 



where j = e l27T / 3 is a cube root of unity, so that its exchange 
properties are very simple. The cyclic exchange (x±, x 2 , x 3 ) — * 
(x 2 , x 3 , Xi) is identical with <fi — > <fr + |7T and the interchange 
x\ <-» x 2 is identical with <fi — > —0. Note that these two permu- 
tations generate the whole symmetric group S3. 
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One easily obtains the interparticle coordinates 

Xjs = y/2 r sin 0, x 23 = y/2 r sin(0 + ?f), (39) 



^31 



y/2r sh 



47T ' 

3 - 



From these relations, one deduces that the value of corresponds 
to the ordering of the three particles as follows, 



O<0<- 

2^ < < 7T 
7T 

7T < < 4- 
o 

7T 7T 

4- < < 5- 
3 3 

5^ < < 27T 



aci > x 2 > x 3 , 

Xl > X 3 > X 2 , 
X 3 > Xi > X 2 , 
X3 > X2 > Xl, 
Xl > X3 > Xl, 
X 2 > Xl > x 3 , 



(40) 



At a bound = n7r/3 with n integer, two particles coincide. On 
the other hand, note the useful identities 

2tt An 1 

sin0 sin(0 + — ) sin(0 + — ) = — - sin 30 (41) 



+ 



1 



1 



9 



sin 



■+4f) sin'(0+4f) sin^30 



(42) 



4.2 Separation of variables 



In the case of N — 3 particles, the wavefunction is completely 
separable as ip — S(A)i?(r)$(0) according to ( [l9|) and (p(i|). 
Using (|^|), the angular eigenvalue equation A$ = A$ reads 



9z/(z/- 1) 
sin 2 30 



(43) 



Since the particles are identical, the angular equation can be 
solved in the irreducible representations of S3, 



S 3 = cm + 2EP 



(44) 



The Bose and Fermi representations are of degree 1, and the two 
mixed representations are equivalent and each of degree 2. Let 
me choose to simultaneously diagonalize the hamiltonian and 
the cyclic permutations. The eigenvalues of a cyclic permutation 
are l,j,j because the cube of a cyclic permutation is merely the 
identity permutation. 

The bosonic (fermionic) states are invariant under the cyclic 
exchanges and they are even (odd) under the interchange of two 
particles, 

cm: $(0 + y) = 4>(0) = $(-0), (45) 

§ : $(0+|^)=$(0)=-$(-0) ! (46) 
The remaining states belong to the mixed representations, 

r$(0), r=j,j 2 . (47) 



F: 4>(0+|) 



The states with r = j and the ones with r = j 2 are in correspon- 
dence under the complex conjugation, and also under the inter- 
change of two particles (doing $(0) — * <&(— 0) and — * —0— |7r) 
in conformity with the degree 2 of the mixed representations. 
Their energy spectrums are identical. 



4.3 Angular quantization 



The differential equation (|43j) may be transformed into the hy- 
pergeometric equation by the redefinition $ = |sin30|" <& and 
the substitution z — cos 2 30. It follows that the angular eigen- 
value equation only has two linearly independent solutions, 



$ x (0) = | s in30|"F( 



cos 2 30), (48) 



$ 2 (0) = I sin 30|" cos 30 F( 



2 ' 2 ; 2' C0S ^' 

(49) 

where fj, — |vA. According to the theory of differential equa- 
tions, only the singularities of the differential equation could 
be singularities of its solutions. These singularities occur at 
= mr/3, when two particles coincide, and they are of alge- 
braic type, namely]^], 



$r(0) = 7 11 |sin30|"F( 
+ 7l2 |sin30| 1 - 1 ' 



v „ f H + v v-fj, l + 2v . 



2 ' 2 



sin 30) 



(50) 



l-/Z-2^1+/i-^3-2f . 2 

— o — . — n — ; — o — ; sm 3 W' 



$2(0) 



721 1 sin 30|" sign (cos 30) 

„,u-\-v v— u, \-\- 2v n 

XF{ 2 ' 2 ;Shl 30) 

1-1/ 



(51) 



+722 1 sin 30| 



sign (cos 30) 
1 — [i — v 1 + n — v 3 — 2v . 2 



sin^ 30), 



with 



7n = 



721 = 



r(i)r(i^) 



r( 



)r(to) 



712 = 



r(|)r(^) 
r(^±^)r(^f^) 



722 



2 I 

r(|)r(^i) 



r( 



)r( 



1 — U+v 



where V is the gamma function. 

At an ordinary point, the eigenfunction is a combination of 
these two independent solutions, 

$(0)=a n $ 1 (0)+^ n $ 2 (0), 0e]ri|,(n+l)|[. (52) 

At a singular point = n7r/3 two particles coincide so that the 
asymptotic behaviour of the wave function must be 



$(0) = I sin30|"/(0) + sin 30 | sin30| _ "g(0) 



(53) 



with / and g continuous, in conformity with (|12|), ( J39| ) and (41). 
The general form of the eigenfunction (|5^) reproduces the alge- 
braic singularities | sin30|" and sin 30 | sin30|~", and the func- 
tions / and g in factor are continuous at = n7r/3 when one 
has 



7n (-)™72i 

712 (— )"722 



Oi n -l 



711 
-712 



(-)"721 
-(-)"722 



a„ 



(54) 

2 by using the relation F(a, b; c; 2) = r(e^a)T(c— 6) ^( a ' b;a + b — c+ 1; 1 — 



*) + (!- zY- a - b T{ f^ c) F{c - a,c - 6;c - a - 



6 — 1:1- 2), and also 

the relation F(a, fe; c; 2) = (1 — z) c ~ a ~ h F(c — a, c — 6; c; 2) in the case of $2 . 
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Note that the functions / and g are then proportional to hyper- 
geometric series: their successive derivatives are also continuous 
at (f> = nir/3. The recurrence relation (|54|) may be rewritten as 



with the transfer matrix 
1 



= T n 



a n - 

Pn- 



T 

1 n 



711722 — 712721 



(55) 



(56) 



711722 + 712721 
-2(-)" 7ll7l2 



2(-)"72l722 

-711722 — 712721 



For a given (ao,/3o)> this recurrence relation determines 
(ai, fix), . . . , (CU5, P5) and thus the wave function is obtained on 
the period [0, 2ir]. However, only the wave functions with pe- 
riod 2ir are physically acceptable, and such functions are only 
for quantized values of A. 

Let me perform the quantization taking advantage of the clas- 
sification into representations of S3. The periodicity requirement 
is then replaced by <E>(</) + |-7r) = r $>((/)) with r = 1 for bosons 
and fermions, and r — j,j 2 for the mixed representation. This 
may be expressed as 



Po 



Po 



(57) 



and has a solution when det^Ti — r) = 0. This last equation 
is responsible for the quantization of A, it may be expanded as 



det T 2 T 1 - r tr T 2 T 1 + r 2 = 0. 
Now, one has det T n = — 1 and 

tr T 2 T X =4 f 7ll722 + 712721 

V 711722 — 712721 , 

Applying the product formula for the gamma function, 

1 COS flTT 



711722 + 712721 



711722 — 712721 



1 



2v cos vtt ' 
1 



1 - 2v 

Finally, the angular eigenvalue A is discretized according to 

^2 



cos fin 



(1 



4r 



■ cos isir, 



y/X = 3fi>0. 



(58) 
(59) 

(60) 
(61) 

(62) 



The solution is immediate by now. Fig.l displays the discrete 
values of vA for the 3! ground states. One clearly sees the 
mirror symmetry and, at v = 1/2, the crossing associated with 
the degenerate solution (|36|). 

The Bose representation is characterized by r = 1 and fi — 
£ + v where £ is a non-negative integer. One has then 712 = 0, 
Pn = 0, a n — a if £ is even; 722 = 0, a n = 0, P n = Po if 
£ is odd. In each case, the remaining hypergeometric function 
degenerates into a Gegenbauer polynomial C/, hence the Bose 
eigenstate basis 



<f>i(<t)) = |sin30| v C£(cos3<£), 



At v = 1, the free states of Fermi statistics are reproduced but 
with some phase shifts at (f> — mr/3. 

The Fermi representation is characterized by r = 1 and fi = 
£+ 1 — v where £ is a non-negative integer. One has then 711 = 0, 
P n =0, a n = ao if £ is even; 721 = 0, a n — 0, p n = Po if I is odd. 
One recovers the Fermi eigenstate basis obtained by Calogero, 

<f>i(<j)) = sin30| siQ3^| _1 'C r j~ , '(cos3^), yf. 



\i = U + 3 - 3i/. 

(64) 

The fermi eigenstates correspond to the Bose ones under the 
mirror symmetry. 

At least, the mixed representation involves the remaining two 
possibilities r = j and r = j 2 . The eigenvalues are solutions of 
the same equation cos 2 fnr = ^ 
they are two-fold degenerate. One obtains 



rr 3 3 . ( COS V1T 

vA = 3£+-=F — arcsm I — - — 

2 7T \ 2 



(65) 



where £ is a non-negative integer. These eigenvalues reproduce 
the non- linear levels in Fig.l, the sign — is responsible for the 
increasing levels and the sign + is for the corresponding ones 
under the mirror symmetry. The eigenstates directly follow from 
"55b and @. 



4.4 On the Bethe ansatz 

In the limit where the harmonic attraction vanishes, the radial 



eigenfunction (20) becomes a Bessel function J^(kr) and the 
relative energy is then fc 2 /2 where k is a positive real. In com- 
parison with a system of free particles, the asymptotic behaviour 
of the eigenfunction 



irkr V 2 4 



(66) 



sets up a phase shift proportional to the ^-dependence of y/X. 
In the case of Bose and Fermi statistics, this phase shift can be 
viewed as a sum of 2-body phase shifts in accordance with the 
Bethe hypothesis ||. In the case of mixed statistics, however, 
it is not the case and the Bethe hypothesis j2Cj is no longer 
relevant. 

5 A numerical evaluation for four 
bodies 

Leaving out the center-of-mass coordinate, I use the Jacobi co- 
ordinates 

1 



(63) 



y = ^(xi-x 2 ), 

x = —=(x! + x 2 - 2x 3 ), 
V6 
1 

(xi + x 2 + x 3 — 3^4). 



c.p. 



(67) 



The relative hamiltonian reads 



~ TT I d v d x - V2d z 



y x - V2z 



(68) 
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where c.p. stands for the cyclic permutations of xi,x 2l x 3 . Ex- 
pressing x, y, z in terms of the spherical coordinates r, 9, cj>, the 
wavefunction may be separated as -0 = E(X)R(r)Y(0,(j)). The 
center-of-mass eigenstates are (|l9| ) and the radial eigenstates 
R = r 3 "R are ©. 

The first angular eigenstates may be evaluated numerically by 
direct diagonalization of the matrix of the angular part of the 
hamiltonian in a basis of spherical harmonics Yi, m . For simplifi- 
cation, the analysis is restricted to the symmetric wavefunctions 
under the exchanges of the particle coordinates x\,x%,X3, that 
is, Yi i7n have to be symmetrized under <fi —> <f>+ |7r and 4> ~* ~4>- 
In such a situation, there are only 4!/3! = 4 ground states. Since 
the potential (x — y/2z)~ 1 (d x — \/2d z ) may be transformed into 
z d z by means of a simple rotation, one has to calculate the 
matrix elements of x d x and z d z . In each case, one has to 
effect a simple Cauchy integration over </> and then the integra- 
tion of a certain polynomial by means of the Gauss-Legendre 

quadrature. Fig. 2 displays the discretized values of JX + i ob- 



tained by the numerical diagonalization of a 222 x 222 matrix 
for different values of v. One distinguishes the 4 ground states 
and some excitements. The Bose groundstate displays a linear 
energy of slope 6. The other energies are not linear but their 
slope is always 2 at v = and v = 1. 



6 Statistical mechanics in the ther- 
modynamic limit 

In this section, the harmonic well has to be understood as a 
long distance regulator alternative to the box For instance, 
expressing the nth cluster coefficient b n = Z n — Z n -\Z\+. . . \Z™ 
in terms of the TV-body partition functions Zn fll7|| , the to — > 
limit can be identified to the infinite box thermodynamic limit 
if the divergent factor n _1 / 2 (/3w) _1 surviving in b n is identified 
to an additive factor, namely VX^, 1 where V is the box volume 
and At = /mkBT is the thermal wavelength. 

In the case of Bose and Fermi statistics, the thermodynamic 
limit of the thermodynamical potential has been obtained by 
various means from the many-body spectrum in a harmonic 
well 



21 



, |23|. The result is in accordance with the ther- 
modynamic Bethe ansatz and thus provides a confirmation of 
the Bethe hypothesis. Interesting enough, the particles obey 
Haldane's generalization of the Pauli principle |3j]. Note that a 
correlation function has recently been analysed |18[ p"9[ . 

In the case of Boltzmann statistics where no exchange sym- 
metry is imposed to the wavefunction, I can calculate the cluster 
coefficients h = Z x , b 2 = Z 2 - \Z\ and b 3 = Z 3 - Z 2 Z X + \Z\ 
from the Boltzmann partition functions Zn = 4n tr N . The 



thermodynamic limit gives b\ 



v\-\ 



AT! 

and the new result 



63 = 



V 1 



3 . COS WK 

arcsm 



(- 

V 7T 



(69) 



This last result is rather different from that of a system of par- 
ticles obeying Haldane statistics || . 



On the multi-species Calogero 
model 



One considers the iV-particle hamiltonian 



H 



N 

E 

2=1 



1 2 2\ 

TUJ X; 



(70) 



where the masses have been eliminated by a suitable redefini- 
tion on the coordinates and on the parameters. The discussion 
displayed between (Jll|) and (Elf) is immediately generalizable to 



the multi-species model if the change v — > i/y is done. However 
the wavefunction is no longer in an irreducible representation of 
Sn- Moreover, the non-unitary transformation 



i,j< 



introduces three-body terms, 



N 



i=l 



i,i< 



(71) 



(72) 

By symmetrizing under the cyclic exchanges of k and reduc- 
ing to the same denominator the term of the sum, one verifies 
that the three-body potential vanishes when Vij — v. 

The possibility to generalize the algebra of Raising and lower- 
ing operators to the multi-species Calogero model by the change 
v — > vij in (^2|) have been studied in [|l2|. It appears that this 
new algebra of operators does not provided physical excitements 
except for center-of-mass and radial excitements. Nevertheless, 
every excitement provided is good at first perturbative order in 

Let me exemplify the multi-species Calogero model by a nu- 
merical and perturbative analysis of the 3-body problem. Only 
the angular part of the eigenvalue problem have to be considered 
again, 



A 



-dl 



V\2{vi2~\) V 23 (V23~1) , ^3l(^31-l) 



< + ¥) 



(73) 



The non-unitary transformation 





V\2 


■ fj. 27r N 




At: 

sin(0+ — ) 


sin <j) 




sin(^+ —) 













H4>) 

gives an operator 

A = -dl + [yxi + v 23 + v 31 ) 2 



27T, 



-2vi2 cot 09^ - 2v 23 cot (4> + — ) 84, 



-2v 3 i cot 



47T 

:*+ T : 

^23^31 



(75) 



^31^12 



sm((j)+ ^-)sm((t) + m sm(4>+ ^)sin0 



3 

^12^23 



sin0sin(0 + ^) 



the matrix elements of which are well defined with the princi- 
pal value regularization. Indeed, the angular wave function is 



7 



periodic of period 2ir and thus it may be expanded in Fourier 
series, 



cte**, 



(76) 



and the matrix elements of A are expressed in the Fourier space 
as simple Cauchy integrals. The final result reads 



(77) 



+fact(£ - £') h-v cvc 

X \Vy&Vz\ + ^31^12 f 



V\2 ^23 J 



■M-1t' 



with sign (0) = and fact (I - (!) = 0, -2, 2 respectively if \£ - 
C\ =0,2,4 modulo 6. Diagonalizing a 361 x 361 matrix, I have 
obtained a numerical evaluation of the first eigenvalues A. Note 
that the energy spectrum is built from the square root of A 
according to (pl|). 

In the case — v, the numerical results reproduce the exact 
spectrum of Fig.l as it should. 

In the case vy^ — and z/ 2 3 = ^31 = v, the first numerical 
values of \f\ are displayed in Fig. 3. The spectrum is nearly 
linear in the quantum number £. However, the computation at 
second perturbative order around v = bears out the presence 
of some irregularities. Since the problem is invariant under 11 <-> 
X2, the non-degenerate perturbative theory may separately be 
used for the unperturbed bases {cos(£(p)} and {sm(£cj))}. For 
instance, starting from the unperturbed eigenstate cos(£4>) with 
I = 1 modulo 3, one finds 



V\ = £-v + { 



* l \ 2 



0(v 3 ) 



(78) 



Since these energies are not linear in the quantum number, there 
is no algebra of Raising and lowering operators in the multi- 
species Calogero model. 

8 Conclusion 

The nature of the wavefunction of the Calogero model has been 
discussed at length. 

As regards the case of Boltzmann statistics, we have revisited 
the Raising and lowering operators and we have obtained some 
new exact eigenf unctions. One would like to derive the complete 
solution of this many-body problem, if possible analytically. 

In the multi-species case, we have shown the impossibility 
of an algebra of Raising and lowering operators. However, the 
question of the integrability of this model remains open. 

At least, it would be desirable to have a better understand- 
ing of the commutation rules between the Raising and lowering 
operators of the one-species model. 
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